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HIGHLIGHTS 


GRAPHICAL  ABSTRACT 


•  p(02)  was  visualized  at  the  GDL  and 
the  upper  channel  wall  of  a  running 
PEFC. 

•  Two  different  oxygen-sensitive 
luminescent  dye  films  were  used. 

•  Different  distributions  of  p(02)  were 
seen  at  the  two  locations  inside  a 
PEFC. 

•  A  numerical  calculation  was  carried 
out  to  understand  the  reactions  in¬ 
side  the  MEA. 

•  Water  distribution  influenced  the 
current  density  and  p(02)  in  the  MEA. 
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Visualization  of  the  oxygen  partial  pressures  was  carried  out  at  the  surface  of  a  gas  diffusion  layer  (GDL) 
for  the  first  time  together  with  the  upper  part  of  the  gas-flow  channel  of  the  cathode  of  a  running 
polymer  electrolyte  fuel  cell  (PEFC)  using  two  different  oxygen-sensitive  luminescent  dye  films.  The 
visualized  distributions  of  the  oxygen  partial  pressures  at  the  GDL  and  the  upper  gas-flow  channel  during 
the  PEFC  operation  were  very  different  in  a  conventional  test  cell.  The  change  in  the  distribution  of  the 
oxygen  partial  pressures  was  observed  by  changing  the  oxygen  utilization,  which  should  be  connected 
with  the  reactive  locations  in  the  membrane-electrode  assembly  (MEA).  A  numerical  calculation  was 
carried  out  to  understand  the  distributions  of  water  and  current  density  inside  the  MEA.  The  water 
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distribution  inside  the  MEA  was  confirmed  to  strongly  affect  the  distributions  of  the  current  density  and 
the  oxygen  /partial  pressure. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Polymer  electrolyte  fuel  cell  (PEFC)  is  a  new  energy  source  ex¬ 
pected  to  be  widely  used  with  its  high  power  density  and  cleanli¬ 
ness.  Although  PEFCs  have  many  advantages  compared  with 
current  internal-combustion  engines,  the  cost  reduction  and  the 
improvement  of  durability/reliability  are  still  needed.  The  devel¬ 
opment  of  new  materials  for  the  cell  components  and  suitable  cell 
designs  are  required.  Adequate  operation  modes  also  need  to  be 
researched.  For  the  improvement  of  PEFCs  via  the  issues  pointed 
out  above,  it  is  important  to  reveal  the  distribution  of  chemical 
species  (O2,  H2,  H2O,  CO2,  H2O2,  etc.)  and  physical  parameters 
(temperature,  current  density,  etc.)  in  a  PEFC  under  the  operating 
conditions  and  to  feedback  the  so-obtained  data  into  the  designs  of 
the  membrane  electrode  assembly  (MEA)  and  the  PEFC  itself.  The 
visualization  inside  the  PEFC,  therefore,  is  essential. 

Recently,  visualization  studies  inside  PEFCs  have  been  focused 
mainly  on  liquid  water  to  understand  the  flooding  and  plugging, 
which  could  occur  in  the  catalyst  layer,  gas  diffusion  layer  (GDL), 
and  gas-flow  channels  to  block  supply  gases.  The  water  visuali¬ 
zations  have  been  carried  out  by  using  various  technologies:  direct 
optical  visualization  1-6  ,  neutron  radiography  imaging  [7-14], 
X-ray  imaging  [15-18  ,  and  magnetic  resonance  imaging  (MRI) 
[19-21].  Those  studies  revealed  that  the  generated  water  plays  a 
key  role  in  the  performance  and  stability  of  the  PEFCs.  Now,  the 
visualization  of  other  chemical  species,  such  as  oxygen,  is 
expected. 

Our  research  group  has  developed  a  nondestructive,  real-time/ 
space  visualization  system  for  oxygen  partial  pressure  (p(02))  at 
the  cathode  of  PEFCs  5,6,22-24]  and  a  direct  methanol  fuel  cell 
(DMFC)  25]  under  the  operating  conditions  using  a  chemical  probe 
for  p(02),  “platinum-tetrakis(pentafluorophenyl)porphyrine 
(PtPP)”  (Fig.  1(a))  [26  .  In  the  previous  papers,  p(02)  was  visualized 
only  on  the  upper  channel  walls.  In  this  paper,  the  distributions  of 
p(02)  at  the  GDL  surfaces,  together  with  that  on  the  upper  channel 
walls,  at  the  cathode  were  visualized  for  the  first  time  using  two 
different  dyes:  regularly-used  PtPP  and  a  new  probe  “platinum- 
tetrakis(pentafluorophenyl)porpholactone  (PtPL)”  (Fig.  1(b)) 
[27,28].  Experimentally,  by  using  those  dyes,  different  p(02)  dis¬ 
tributions  were  visualized  between  at  the  GDL  surface  and  at  the 
upper  position  of  a  gas-flow  channel. 
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Fig.  1.  (a)  Chemical  structure  of  platinum-tetrakis(pentafluorophenyl)porphyrine 
(PtPP).  (b)  Chemical  structure  of  platinum-tetrakis(pentafluorophenyl)porpholactone 
(PtPL).  (c)  Absorption  spectrum  and  luminescence  spectra  of  PtPP  and  PtPL. 


For  understanding  the  oxygen  reduction  reaction  (ORR)  in  the 
MEA,  or  the  reaction  current  distribution,  we  also  carried  out  nu¬ 
merical  and  simulation  inside  the  PEFC.  To  evaluate  the  perfor¬ 
mance  of  PEFCs  speedily  and  accurately,  our  research  group  has 
developed  a  simulation  tool  including  the  macroscopic  parameters 
of  gas  and  heat  transport  in  flow  field  and  MEA  characteristics  such 
as  electrochemical  kinetics  in  the  catalyst  layer  and  proton/water 
transport  in  the  electrolyte  membrane  [29  .  In  this  paper,  experi¬ 
mental  and  simulation  data  are  carefully  compared  to  demonstrate 
a  new  evaluation  strategy  for  understanding  the  reaction  distri¬ 
butions  inside  a  PEFC.  The  simulation  results  gave  clear  explana¬ 
tions  for  the  change  of  the  p(02)  distributions  under  different  cell- 
operation  conditions. 

2.  Experimental 

2  A.  Oxygen  sensitive  dye  films  and  see-through  PEFC 

Toluene  solutions  of  PtPP  (Fig.  1(a))  and  PtPL  (Fig.  1(b))  were 
mixed  with  an  oxygen-permeable  polymer,  poly(l-trimethylsilyl- 
1-propyne)  (PMSP)  at  room  temperature  as  previously  reported 
[27,28].  Absorption  and  emission  spectra  of  PtPP  and  PtPL  are 
shown  in  Fig.  1(c).  Both  PtPP  and  PtPL  have  the  absorption  peaks 
approximately  at  400  nm,  but  have  different  emission  peaks  at  650 
and  750  nm,  respectively.  Therefore,  p(02)  at  two  positions  inside  a 
cell  can  be  visualized  by  using  a  single  light  source  (407  nm). 

Commercial  MEAs  with  the  Pt  loadings  of  0.4  mg  cm-2  both  on 
the  anode  and  the  cathode  were  used.  The  catalyst  area  was 
52  mm  x  52  mm.  For  the  electrolyte  membrane  with  a  thickness  of 
30  pm,  perfluorosulfonic  acid  ionomers  were  employed.  On  a  gas 
diffusion  layer  with  a  thickness  of  approximately  200  pm,  a 
microporous  layer  (MPL)  with  a  thickness  of  several  micrometers 
was  attached. 

A  see-through  cell  with  a  transparent  window  was  used  for 
visualizing  p(02)  in  the  PEFC  under  the  operating  condition.  The 
endplate  of  a  “Standard  Cell  of  Japan  Automobile  Research  Institute 
(JARI  standard  cell)”  [30],  a  commonly-used  cell  employed  for  the 
development  of  PEFCs,  at  the  cathode  side  was  replaced  by  an  acrylic 
plate  to  visualize  inside  the  cell  (Fig.  2(a))  [24,31  .  PtPL  was  sprayed 
on  a  perfluoroalkoxy  (PFA)  film  (film  thickness  =  2  pm),  and  the  film 
was  inserted  between  the  endplate  and  the  current  collecting  ribs.  In 
this  way,  p(02)  on  the  upper  wall  of  the  gas-flow  channel  was  visu¬ 
alized.  On  the  other  hand,  PtPP  was  sprayed  uniformly  on  the  GDL 
surface  at  the  cathode  with  the  average  thickness  of  4  pm.  The 
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Fig.  2.  (a)  Photograph  of  a  see-through  cell,  (b)  Schematic  drawing  of  a  flow  field  with 
two  dye  films. 
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locations  of  the  dye  films  are  schematically  illustrated  in  Fig.  2(b).  The 
channel  depth  and  width  of  the  cell  were  both  1  mm,  therefore,  the 
distance  between  the  PtPP  and  PtPL  films  was  also  1  mm. 

2.2.  Optical  settings  for  the  visualization  of  p(02) 

The  optical  setting  for  visualization  of  p( O2)  was  essentially  the 
same  as  that  reported  in  previous  papers  [5,22,24,25  .  Fig.  3  shows 
a  schematic  drawing  of  the  optical  system.  A  diode  laser  was  used 
as  a  407-nm  light  source,  and  the  excitation  light  was  diffused  by 
diffuser  plates  and  distributed  uniformly  onto  the  see-through  cell. 
The  emission  lights  from  PtPP  at  the  GDL  surface  and  PtPL  at  the 
upper  channel  wall  passed  through  an  acrylic  endplate  and  were 
separated  by  a  hot  mirror.  Then,  the  lights  were  separately  captured 
with  two  CCD  cameras  (500  x  500  pixels,  1  pixel  =  100  x  100  pm) 
in  a  dark  room.  CCD  camera  1  in  Fig.  3  with  a  band-pass  filter 
(620  <  A  <  700  nm)  in  front  was  used  for  detecting  P(02)  on  the  GDL 
surface,  whereas  CCD  camera  2  with  a  cutting  filter  (A  >  720  nm)  for 
detecting  p(02)  on  the  upper  channel  wall.  In  this  way,  the  images 
on  the  GDL  and  on  the  upper  wall  of  the  gas-flow  channel  were 
separately  and  simultaneously  obtained. 

2.3.  Visualization  of  p(0 2) 

Prior  to  the  power  generation,  the  calibration  curves,  known  as 
Stern-Volmer  plots,  of  PtPP  and  PtPL  were  obtained  by  the  same  way 
as  reported  previously  [22,24,25]  under  the  mixed  gases  of  oxygen 
and  nitrogen  (oxygen  concentration  =  0-25%)  at  80  °C  supplied  to 
both  cathode  and  anode  gas-flow  channels  at  the  gas-flow  rate  of 
50  mL  min-1.  The  cell  operation  was  performed  at  80  °C,  with  the 
humidified  air  and  hydrogen  (60%  RH)  being  supplied  to  the  cathode 
and  the  anode,  respectively,  at  300  mL  min-1.  Oxygen  uti- 
lization(Uo2)  was  set  at  0,  20,  30,  60,  and  80%.  During  the  operation, 
the  laser  light  was  irradiated  onto  the  cell,  and  the  light  emission 
from  the  GDL  surface  and  the  upper  channel  wall  were  separately 
captured  by  CCD  cameras.  The  time  and  space  resolutions  were 
250  ms  and  100  pm,  respectively.  Under  the  steady  state  condition, 
32  emission  images  were  continuously  captured  to  be  averaged  on 
the  computer.  Distributions  of  the  emission  were  converted  to  the 
p( O2)  images  on  the  upper  channel  wall  and  on  the  GDL  surface  using 
the  Stern-Volmer  plots  acquired  previously  for  both  dyes. 

3.  Numerical  modeling 


MEA,  are  nonlinearly  coupled  [32  .  These  phenomena  require  a 
large  amount  of  efforts  to  solve  many  equations  by  fully  three- 
dimensional  fundamental  equations  based  on  Computational 
Fluid  Dynamics  (CFD),  so  it  is  necessary  to  apply  simplified  models 
and  cut  the  calculation  cost  down.  In  this  study,  macroscopic 
models  are  applied  for  the  phenomena  inside  the  MEA,  such  as 
electrochemical  reactions  and  proton/water  transport.  Further¬ 
more,  these  models  are  coupled  with  multi-dimensional  heat  and 
mass-transport  equations,  in  which  effective  parameters  such  as 
gas  diffusivity  and  permeability  in  the  GDL/MPL  and  the  equivalent 
diameter  of  the  flow  channel  are  included.  These  parameters  are 
emphasized  with  superscript  “eff’  in  the  following  equations. 

3.2.  Models  for  mass  and  heat  transport 

In  this  study,  the  mass  and  heat  transports  in  the  flow  field  are 
modeled.  Reynolds  number  in  the  cathode  channel  under  the  pre¬ 
sent  experimental  conditions  was  first  calculated  to  confirm  the 
laminar  flow  in  the  channel.  In  the  flow  channel,  the  one¬ 
dimensional  gas,  liquid  water  and  heat  transport  equations,  in 
which  the  pressure  drop  is  assumed  to  be  of  the  laminar  flow,  are 
applied.  In  the  GDL  with  an  MPL,  the  fully  three-dimensional  gas, 
liquid  water,  and  heat  transport  equations  with  multi-component 
gas  diffusion  are  applied,  in  which  the  pressure  drop  is  assumed  to 
follow  the  Darcy's  law.  The  advantage  of  these  simplified  models  for 
the  flow  calculation  is  a  substantial  reduction  of  calculation  cost  in 
comparison  with  CFD.  Therefore,  these  models  can  simulate  a  large- 
scale  cell  performance  in  a  short  period  (about  half  a  day  using  a 
most  recent  personal  computer).  These  equations  are  shown  below. 

•  conservation  equation  of  mass: 


[pg(  1  -  s)]  +  V-  [pg(  1  -  s)ug]  -  rhpS  +  mjr*2° 

*  (1) 

k  + v-(p/su/)  =  -™f2° 

•  equation  of  pressure  loss  in  the  flow  field  (based  on  Darcy's 
law): 


Inside  PEFCs,  multi-physical  and  multi-scale  phenomena,  such 
as  gas/heat  transports  in  the  flow  channel  and  the  GDL  with  an 
MPL,  electrochemical  reactions,  and  proton/water  transports  in  the 
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In  the  flow  channel,  the  permeability,  I<,  is  replaced  by  the 
parameter  with  equivalent  diameter,  de ,  of  a  laminar  pipe  flow,  and 
the  capillary  pressure,  pc,  as  a  function  of  water  saturation  is 
assumed  to  be  zero. 


(3) 


•  conservation  equation  for  chemical  species: 
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Fig.  3.  Schematic  representation  of  the  optical  settings  for  the  p02visuaFzation. 
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where  m^°  is  the  mass  production  by  phase-change  that  is 
formulated  as  follows: 


Heat  production  by  the  over-potential  is  calculated  by  the 
following  equation. 


m 


h2o 

gi 


(5) 


•  conservation  equation  of  energy 


(1  -  r)  —  pshs  +  £—  [pg{  1  -  s)hg  +  Pishi] 

+  V-  [Pg(  1  -  s)hgUg  +  Pishiu,}  =  V-  (>ffV7)  +  <j(m) 


3.2.  Models  for  electrochemical  reactions  and  for  transports  of 
protons  and  water 


Electrochemical  reactions  and  water  transports  in  the  MEA  are 
treated  with  parametric  models,  which  involve  the  partial  pressure 
of  reactants  and  the  temperature  at  the  closest  mesh  point  to  the 
MEA  as  representative  parameters.  Electrochemical  reactions  are 
assumed  to  follow  the  Butler- Volmer  equation  and  occur  at  the 
interface  between  the  MPL  and  the  electrolyte  membrane.  The 
catalyst  layer  is  treated  as  a  planar  surface. 


*  —  l0(/<) 


exp 
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Here,  0^  is  catalyst  utilization  which  depends  on  the  MEA 
characteristics.  Ohmic  over-potential  due  to  a  proton  transport  in 
the  electrolyte  membrane  is  described  by  the  following  equations 
according  to  the  model  proposed  by  Nguyen  and  White  [33  . 


(9) 


(10) 


Each  mass  flux  of  the  chemical  species  generated  by  the  elec¬ 
trochemical  reaction  is  given  by  the  following  expressions.  As  for 
the  water  transport  through  the  electrolyte  membrane,  both  the 
electro-osmotic  drag  and  the  water  diffusion  are  considered. 


•  Mass  flux  by  electrochemical  reactions 


Q(m) 


i 

f(m) 


r!a(a)  T-  Va(c )  +  Vr 


(13) 


For  solving  the  continuum  Equations  1-13  under  the  bound¬ 
ary  conditions  shown  in  Table  1,  it  is  necessary  to  determine  the 
macroscopic  parameters.  These  parameters  could  be  determined 
directly  using  the  results  of  a  micro-scale  simulation,  for  example, 
molecular  dynamics  simulation  in  the  catalyst  layer  and  the 
electrolyte  membrane  and  fully  Navier-Stokes  simulation  in  the 
flow  field.  In  this  study,  however,  they  are  adjusted  to  reproduce 
quantitatively  the  measured  polarization  curves  of  a  small-sized 
cell  and  pressure  drop  in  the  flow  channel  under  several  condi¬ 
tions  for  different  relative  humidities,  oxygen  gas  concentrations 
and  cell  temperatures.  The  reactant  gas  utilization  is  very  low 
(less  than  1%)  in  order  to  obtain  the  sufficient  liquid  water 
removal  in  the  flow  field  and  the  uniform  current  distribution 
[29].  The  final  fitting  parameters  using  in  this  study  are  shown  in 
Table  2. 


4.  Results  and  discussion 

Fig.  4(a)  and  (b)  show  the  Stern-Volmer  plots  of  PtPP  and  PtPL 
dye  films,  respectively,  placed  in  a  see-through  cell.  pref  is  the  ox¬ 
ygen  partial  pressure  and  7ref  the  emission  intensity  obtained  in  air. 
Both  the  Stern-Volmer  plots  of  PtPP  and  PtPL  appeared  nearly 
linear.  The  plots  were  fit  with  cubic  polynomials,  which  were  used 
for  the  calibration  curves.  The  performance  of  the  see-through  cell 
with  PtPL  and  PtPP  dyes  coated  inside  was  examined  before  the 
visualization.  Fig.  5  shows  the  current  density-voltage  (/-V)  curves 
obtained  in  a  see-through  cell  with/without  dye  films.  The  I-V 
curve  of  a  see-through  cell  without  dyes  (blue  line)  was  approxi¬ 
mately  identical  to  that  of  a  JARI  standard  cell.  With  a  PtPL  dye  film 
on  the  upper  channel  wall  (red  line),  the  I-V  curve  overlapped  with 
that  without  dyes.  Additional  coating  of  a  PtPP  dye  film  on  the  GDL 
surface  slightly  lowered  the  voltage  at  a  high  current  density 
above  0.7  V,  probably  because  of  the  smaller  gas  diffusivity  at  the 
GDL  surface.  Overall,  the  cell  performance  was  high  even  with  dye 
films. 

Fig.  6  show  the  time-resolved  visualization  images  of  the  dis¬ 
tributions  of  p(02)  at  the  upper  channel  wall  (Fig.  6(a)-(e))  and  at 
the  GDL  surface  (Fig.  6(f) — (j))  simultaneously  obtained  at  the 
cathode  during  the  cell  operation.  The  color  of  the  images  corre¬ 
sponds  to  p(02).  Air  was  supplied  from  the  upper  left  and  was 
exhausted  from  the  lower  right.  Fig.  6(a)  and  (f)  show  the  steady- 
state  p(o2)  at  Uq2  =  20%.  Then,  U02  was  abruptly  changed  from 
20%  to  80%,  and  the  data  acquisition  was  carried  out  with  a  time 
resolution  of  250  ms.  At  the  upper  channel  wall  at  0.25  s  (Fig.  6(b)), 
in  which  the  signals  were  accumulated  between  0  and  0.25  s,  no 
distinctive  change  was  seen.  However,  on  the  GDL  surface 
(Fig.  6(g)),  the  change  in  p(02)  was  clearly  observed,  showing  the 


• j< 

mp 


k  0 
— W  — 

dZ 


•  Water  transport  flux 


=  -w 


H?0 


0 

dZ 


P(m)neff 
EW  w 


(11) 


(12) 


Table  1 

Boundary  conditions. 


Parameters 

Values 

Anode  inlet 

Hydrogen:  300  mL  min  1 

Relative  humidity:  60%  RH 

Temperature:  353.15  K 

Cathode  inlet 

Air:  300  mL  min  1 

Relative  Humidity:  60%  RH 

Temperature:  353.15  K 

Anode/cathode  outlet 

Pressure:  101.325  kPa  (abs) 

Cell  temperature 

353.15  K  (on  the  upper/lower  boundary  of  a  cell) 
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Table  2 

Model  parameters  for  MEA  and  GDL. 

Parameters  Values 

Butler— Volmer  equation  (electrochemical  io(a),ref  =  1.0  x  10°  A  cur2 

reaction)  a^a)  =  0.25,  =  0.5 

l’o(c),ref  =  1.0  x  10-5  A  cm  2 
«fc)  =  1.35,  7(c)  =  0.9 


Catalyst  utilization 

^(a)  —  L0,  $(c)  —  Qy(c) 

Proton  conductivity  in  the  PEM 

®$)  =  0.5A(m)  -  0.326 

Osmotic-drag  and  water  diffusivity 

Springer  et  al.  [34] 

Porosity 

cgdl  =  0.8,  £mpl  =  0.4 

Permeability 

/<gdl  =  4.0  x  10'12  m2 

/Cmpl=  1.0  x  10~13  m2 

Effective  diffusivity 

DffF=Dr[£(1-s)]25 

Equivalent  diameter  of  flow  channel 

de  =  1.0  x  10  3m 

Phase  change  rate 

wgi  =  1.0  x  102  s'1 

lowering  of  p(02)  near  the  outlet.  Notably,  the  decrease  of  p(02)  was 
not  along  the  channel  length,  but  directly  from  the  inlet  to  the 
outlet,  as  a  contour  line  (dashed  line)  around  11  kPa  shows.  At  0.5  s, 
the  change  in  p(02)  became  larger  on  the  GDL  surface  (Fig.  6(h)),  and 
on  the  upper  channel  wall,  the  change  in  P(02)  also  appeared. 
(Fig.  6(c)).  On  the  upper  flow  channel,  the  influence  of  the  cross 
flow  was  much  smaller  than  that  on  the  GDL;  a  contour  line  was 
rather  parallel  to  the  flow  channel.  After  1.0  s,  the  distributions  of 
p(o2)  were  still  observed  to  be  changing.  It  took  several  minutes 
before  the  p(02)  images  were  apparently  stabilized  both  on  the  GDL 
and  the  upper  channel  wall. 

The  steady-state  distributions  of  p0l  are  presented  in  Fig.  7.  The 
exposure  time  for  each  image  was  8.0  s  Fig.  7(a)-(d)  shows  the 
distributions  of  p(02)  at  the  upper  channel  wall,  whereas 
Fig.  7(e)-(h)  at  the  GDL  surface  at  the  cathode  during  the  cell 
operation.  The  visualizations  were  carried  out  at  U02  =  0,  30,  60 
and  80%.  In  our  previous  papers,  we  have  reported  the  visualiza¬ 
tion  of  liquid  water  in  a  DMFC  [25]  and  on  the  GDL  of  a  PEFC  [6]. 
Even  the  simultaneous  visualizations  of  p(C>2)  and  liquid  water 
have  been  carried  out  [5  .  During  the  visualization  of  p( O2)  in  the 
present  study  in  a  single-serpentine  cell,  no  liquid  water  was 
observed  in  the  cathode  gas-flow  channel. 

In  the  case  of  Uq2  =  0%  (0  A  cm-2),  although  the  cell  was  at  the 
open  circuit  with  no  power  generation,  p(C>2)  at  both  the  GDL 
surface  and  the  upper  gas-flow  channel  was  seen  to  decrease 
uniformly  from  15  to  14  kPa  due  to  the  total  pressure  loss  along 
the  gas-flow  channel.  When  the  cell  was  operated  at  Uq2  =  30% 
(0.22  A  cm-2),  p(02)  uniformly  decreased  by  the  ORR  along  the 
gas-flow  channel  from  the  inlet  to  the  outlet  as  seen  in  Fig.  7(b). 


Current  Density  /  A  cm  2 

Fig.  5.  Current  density-voltage  curves  in  a  see-thorough  cell  without  dyes  (blue  line), 
with  a  PtPL  film  on  the  upper  flow  channel  (red  line),  and  with  a  PtPL  film  on  the  upper 
flow  channel  and  PtPP  on  the  GDL  (green  line).  Cell  temperature  =  80  °C;  relative 
humidity  =  60%  (air  and  H2);  a  constant  gas  flow  rate  for  hydrogen  and 
air  =  300  mL  min'1.  (For  interpretation  of  the  references  to  color  in  this  figure  legend, 
the  reader  is  referred  to  the  web  version  of  this  article.) 

However,  p(C>2)  at  the  GDL  surface  was  different;  p(C>2)  signifi¬ 
cantly  decreased  in  the  middle  of  the  gas-flow  channel  with  a 
change  of  colors  from  red  to  green  as  seen  in  Fig.  7(f).  Similarly, 
atUo2  =  60  and  80%,  the  distribution  of  p( O2)  at  the  upper 
channel  wall  was  uniformly  decreased,  but  the  distributions  of 
p(02)  at  the  GDL  surface  were  not  identical  to  those  at  the  upper 
channel  wall.  The  influence  of  the  air  flow  through  the  GDL  is 
clearly  seen,  which  makes  a  p(C>2)  pattern  on  the  right  side  on  the 
GDL  surface,  for  example,  slightly  higher  than  on  the  left  side, 
notably  observed  in  Fig.  7(f  .  For  emphasizing  the  difference  of 
p(02)  at  the  GDL  surface  and  at  the  upper  channel  wall,  A 
p(02)  =  (p(02)  at  upper  channel  wall)  -  (p(02)  at  GDL)  was 
calculated  from  Fig.  7  and  plotted  in  Fig.  8.  At  Uq2  =  0%,  A  p( O2) 
was  nearly  equal  to  0  kPa  at  any  location  as  expected.  At 
Uq2  =  30%,  A  p( O2)  is  larger  in  the  latter  half  of  the  MEA.  At 
Uq2  =  60%,  A  p( O2)  is  the  largest  in  the  upper  middle  area  of  the 
MEA,  and  at  80%,  A  p(C>2)  is  the  largest  at  the  upper  location  of 
the  MEA.  Therefore,  the  reactive  area  on  the  MEA  is  expected  to 
differ  at  different  reaction  conditions. 

To  examine  the  data  in  more  details,  p(C>2)  measured  at  the 
upper  channel  wall  and  the  GDL  surface  are  plotted  along  the 
channel  length.  When  Uq2  was  0%  at  the  open-circuit  voltage 
(Fig.  9),  the  average  p(02)  value  at  the  upper  channel  wall  (red 
line)  was  the  same  as  that  at  the  GDL  surface  (blue  line).  p(02)  at 
the  different  locations  decreased  gradually  from  the  inlet  (15  kPa) 
to  the  outlet  (14  kPa),  because  of  a  pressure  loss.  Periodic,  sharp 
spikes  were  seen  both  on  the  upper  channel  wall  and  the  GDL  at 
the  turns  of  the  serpentine  channel,  but  the  spikes  were  much 
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Fig.  4.  Stern-Volmer  plots  for  PtPL  (a)  and  PtPP  (b)  in  the  mixtures  of  oxygen  and  nitrogen.  Cell  temperature  =  80  °C.  The  emission  intensity,  /,  was  measured  for  different  oxygen 
partial  pressures,  p(02).  /ref  is  the  emission  intensity  at  pre f,  the  oxygen  partial  pressure  in  air. 
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Fig.  6.  Time-resolved,  simultaneous  visualization  of  the  oxygen  partial  pressures  at  the  upper  channel  wall  (a)-(e)  and  on  the  surface  of  GDL  (f)-(j)  in  an  operating  PEFC.  The 
oxygen  utilization  was  changed  from  20%  to  80%  at  0  s.  Cell  temperature  =  80  °C;  relative  humidity  =  60%  (air  and  H2);  a  constant  gas  flow  rate  for  hydrogen  and 
air  =  300  mL  min-1.  Time  resolution  was  250  ms.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


larger  on  the  GDL.  Those  spikes  at  the  turns  will  be  later  explained 
by  the  numerical  simulation.  At  Uq2  =  30,  60,  and  80%,  p( O2) 
measured  at  the  upper  channel  wall  (Fig.  10(b)-(d))  and  the  GDL 
surface  (Fig.  10(f)-(h))  are  plotted  with  gray  lines  along  the  gas- 
flow  channel  length  from  the  inlet  to  the  outlet  (1400  mm  in 
length).  At  U02  =  30%,  p(02)  at  the  upper  channel  wall  gradually 
decreased  from  the  inlet  (15  kPa)  to  the  outlet  (10  kPa)  (Fig.  10(b)), 
whereas  p(02)  on  the  GDL  surface  was  clearly  lower  than  that  at 
the  upper  channel  wall  (Fig.  10(f))  at  the  latter  half  of  channel 
(500-1400  mm).  Periodic  spikes  were  seen  as  in  the  case  at 
Uo2  =  0%  (Fig.  9).  The  change  in  p(02)  at  the  GDL  in  Fig.  10(f)  is 
large  in  the  middle  of  the  length  of  the  gas  channel.  From  the 
point  approximately  at  500  mm  to  the  outlet,  p(02)  at  the  GDL  is 
always  lower  than  that  at  the  upper  channel  wall  (Fig.  10(f)), 
showing  that  the  oxygen  might  not  be  efficiently  supplied  to  the 
GDL  surface  from  the  channel  and  that  some  part  of  air  in  the 
channel  could  be  simply  blown  out  at  the  upper  part  of  the  gas 
channel  to  the  outlet  due  to  the  depth  of  the  channel,  1  mm,  and 
to  a  relatively  high  flow-rate  of  air  (Fig.  2(b).  At  U02  =  60%,  p(02)  at 


the  upper  channel  wall  was  gradually  decreased  from  the  inlet 
(15  kPa)  to  the  outlet  (5  kPa)  (Fig.  10(c)),  whereas  p(02)  at  the  GDL 
surface  steeply  lowered  around  300-500  mm  due  to  the  increased 
ORR  (Fig.  10(f)).  In  the  case  at  the  highest  U02  of  80%,  p(02)  at  the 
GDL  surface  was  drastically  decreased  near  the  inlet  of  the  gas- 
flow  channel  around  200  mm  (Fig.  10(h)).  On  the  other  hand, 
p(02)  at  the  upper  channel  wall  decreased  smoothly  from  the  inlet 
(15  kPa)  to  the  outlet  (2  kPa)  (Fig.  10(d)).  Comparing  the  results  at 
Uq2  =  30,  60,  and  80%  in  Fig.  10,  A  p(02)  was  confirmed  to  move 
forward  along  the  gas-flow  channel  (Fig.  8).  At  a  low  U02  of  30%, 
the  electrolyte  membrane  takes  the  water  vapor  formed  at  the 
upper  air  stream,  which  could  mainly  be  accumulated  at  the  lower 
part  of  the  gas-flow  channel  via  the  GDL  into  the  membrane  and 
the  catalyst  layer.  The  proton  conductivity  of  the  electrolyte 
membrane  thus  should  be  improved,  resulting  in  an  enhanced 
ORR  rate  in  this  region.  At  a  higher  U02  of  60%,  a  larger  amount  of 
water  vapor  was  generated  than  at  Uq2  =  30%  closer  to  the  inlet, 
then  the  proton  conductivity  and  the  reaction  might  have 
increased  at  the  middle  part  of  the  flow  length.  When  U02  =  80%, 


Fig.  7.  Oxygen  partial  pressures  at  the  upper  channel  wall  (a)-(d)  and  on  the  surface  of  GDL  (e)-(h)  visualized  experimentally  in  an  operating  PEFC.  Cell  temperature  =  80  °C; 
relative  humidity  =  60%  (air  and  H2);  a  constant  gas  flow  rate  for  hydrogen  and  air  =  300  mL  min-1;  oxygen  utilization  =  0  (a)(e),  30  (b)(f),  60  (c)(g),  and  80%  (d)(h),  corresponding 
to  overall  current  densities  at  0,  0.22,  0.43  and  0.58  A  cm  2,  respectively.  32  images  were  accumulated  for  each  set  of  conditions. 
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Fig.  8.  A  p(02)  ((p(02)  at  upper  channel  wall)  -  (p(02)  at  GDL))  calculated  from  Fig.  3.  Cell  temperature  =  80  °C;  relative  humidity  =  60%  (air  and  H2);  a  constant  gas  flow  rate  for 
hydrogen  and  air  =  300  mL  min-1;  oxygen  utilization  =  0  (a),  30  (b),  60  (c),  and  80%  (d).  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to 
the  web  version  of  this  article.) 


liquid  water  should  exist  in  a  wide  range  of  the  catalyst  layer  and 
the  membrane,  and  the  reaction  rate  should  be  the  highest  near 
the  inlet  with  a  decent  amount  of  water.  On  the  other  hand,  p( O2) 
at  the  upper  flow-channel  wall  did  not  show  a  steep  change 
because  of  less  influence  from  the  MEA.  The  difference  in  the  p(02) 
changes  between  these  two  layers  is  attributed  partly  to  the  non- 
optimized  flow-channel  design  and  could  be  reduced  by 
increasing  the  gas  flux,  e.g.,  by  reducing  the  rib  height  of  the  gas- 
flow  channel  even  at  the  same  Uq2  (in  practical  applications,  the 
depth  of  the  channels  is  normally  less  than  0.5  mm).  The  change  of 
the  channel  design  might  result  in  a  higher  cell  performance  due 
to  the  smaller  over-potential  for  the  ORR  by  increasing  the  mass 
transfer. 

In  order  to  understand  those  processes  more  profoundly,  a  nu¬ 
merical  simulation  was  carried  out  using  a  fuel-cell  model 
described  in  Numerical  Modeling  section.  Simulated  p(C>2)  both  at 
the  upper  channel  wall  and  at  the  surface  of  GDL  are  plotted  in  dark 
along  the  gas-flow  channel  from  the  inlet  to  the  outlet  in  Fig.  10  for 
the  comparison  with  the  experimental  data  shown  in  gray.  The 
distributions  of  p(02)  were  quantitatively  in  good  agreement  with 
experimental  data  at  high  U02,  especially  at  80%,  both  on  the  upper 
channel  wall  and  the  GDL  On  the  other  hand,  at  Uq2  =  30%,  the 
distribution  of  p(02)  was  different  from  the  experimental  result, 
especially  at  the  downstream  area  on  the  GDL  (Fig.  10(f)).  These 
results  indicate  that  the  oxygen  consumptions  by  the  ORR  esti¬ 
mated  at  the  numerical  simulation  does  not  correspond  to  the 
actual  values  at  low  U02. 


Fig.  9.  Oxygen  partial  pressures  at  the  open-circuit  voltage  at  the  upper  channel  wall 
(red  line)  and  on  the  GDL  surface  (blue  line)  along  the  gas-flow  channel  length  from 
the  inlet  (0  mm)  to  the  outlet  (1400  mm).  Experimental  data  were  obtained  from  the 
visualization  data  in  Fig.  3.  Cell  temperature  =  80  °C;  relative  humidity  =  60%  (air  and 
H2);  a  constant  gas  flow  rate  for  hydrogen  and  air  =  300  mL  min-1.  (For  interpretation 
of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version 
of  this  article.) 


The  distributions  p(02)  on  the  cathode  catalyst  layer,  the  local 
current  density  on  the  MEA,  and  the  water  vapor  activity  in  the 
polymer  electrolyte  membrane  derived  from  the  numerical  simu¬ 
lations  are  shown  in  Fig.  11(a) — (i),  respectively,  at  different  Uq2, 
which  can  be  compared  to  the  experimental  results  shown  in  Fig.  7. 
For  p( O2)  (Fig.  ll(a)-(c))  and  the  water  vapor  activity 
(Fig.  11(g) — (i)),  the  absolute  partial  pressures  and  the  relative  hu¬ 
midities  are  colored,  respectively,  whereas  for  the  local  current 
density,  the  color  shows  the  difference  from  the  average  values  in 
percent.  At  Uq2  =  30%,  p(02)  (Fig.  11(a))  decreased,  and  the  current 
density  (Fig.  11(d))  and  the  water  activity  (Fig.  11(g))  increased 
smoothly  along  the  gas-flow  channel  from  the  inlet  to  the  outlet. 
However,  the  p(02)  values  experimentally  visualized  on  the  GDL 
(Fig.  7(f))  decreased  significantly  in  the  middle  of  the  gas-flow 
channel  (Fig.  10(f))  and  deviated  from  those  simulated  (Fig.  11(a)), 
because  of  the  large  oxygen  consumption.  This  discrepancy  and 
insufficiency  of  the  numerical  model  are  probably  caused  by  the 
insufficient  estimation  of  the  water  vapor  transport  through  the 
MEA,  which  directly  affects  the  proton  conductivity  in  the  catalyst 
layer  and  the  electrolyte  membrane.  According  to  the  experimental 
results,  the  actual  hydration  inside  the  MEA  under  low  relative 
humidity  and  low  current  density  should  be  larger  than  that  esti¬ 
mated  by  the  simulation.  Recently,  we  reported  that  the  water 
distributions  across  the  electrolyte  membranes  during  the  cell 
operation  are  very  complex  [35,36].  Therefore,  the  experimental 
data  of  the  water  distribution  inside  the  membrane  are  essential  to 
improve  the  numerical  simulation.  At  U02  =  60  and  80%,  the  MEA 
should  be  much  hydrated  by  produced  liquid  water  (Fig.  11(h)  and 
(i))  during  the  power  generation  (Fig.  11(e)  and  (f)),  and  the  proton 
conductivity  accordingly  increased.  Therefore,  the  effect  of  water 
vapor  transport  could  be  practically  ignored.  At  higher  U02,  the 
distributions  of  water  activity  in  the  MEA  and  thus  p(02)  data  were 
better  simulated  (Fig.  10(g)  and  (h)),  and  the  peaks  of  current 
density  were  numerically  reproduced  to  shift  closer  to  the  inlet  as 
Uq2  increased  (Fig.  ll(e)-(f)).  The  oxygen  consumption  model  at 
different  Uq2  suggested  by  the  experimental  data  are  computa¬ 
tionally  supported. 

In  the  experimental  data  of  p(C>2)  in  Figs.  8  and  10,  especially  in 
the  data  on  the  GDL,  sharp  spikes  are  seen  both  upwards  and 
downwards.  The  positions  of  the  spikes  were  well  simulated;  the 
spikes  are  observed  at  turns  of  a  serpentine  gas-flow  channel.  The 
sudden  decrease  and  increase  of  p(02)  can  be  explained  by  the  air 
flow  through  the  GDL,  short-cutting  the  flow  channel.  The  intensity 
of  the  spikes  is  as  large  as  several  kilopascals,  which  should  be  large 
enough  to  have  an  influence  on  the  cell  performance. 

In  this  paper,  we  visualized  p(02)  on  the  GDL  surface  for  the  first 
time,  and  compared  the  data  with  those  on  the  upper  flow  channel. 
By  using  our  visualization  method,  we  demonstrated  a  capability  of 
providing  better  operating  conditions  during  the  power  genera¬ 
tions  and  keys  to  a  design  of  better  MEAs  and  flow  fields  of  fuel 
cells.  Numerical  simulation  was  carried  out  for  the  visualization 
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Fig.  10.  Oxygen  partial  pressures  at  the  upper  channel  wall  (b)-(d)  and  on  the  surface  of  the  GDL  (f)-(h)  along  the  gas-flow  channel  length  from  the  inlet  (0  mm)  to  the  outlet 
(1400  mm).  Experimental  and  numerical  calculation  data  are  shown  in  gray  and  dark,  respectively.  Experimental  data  were  obtained  from  the  visualization  data  in  Fig.  3.  Oxygen 
utilization  =  30  (b)(f),  60  (c)(g),  and  80%  (d)(h). 


system.  The  comparison  of  experimental  and  simulation  data  was 
also  confirmed  to  be  useful  in  understanding  the  distributions  of 
various  physical  and  chemical  properties  inside  the  fuel  cell.  In 
order  to  improve  the  accuracy  of  the  simulation,  the  water  trans¬ 
port  model  used  here  must  be  improved.  For  this  purpose,  under¬ 
standing  the  water  distribution  inside  the  membrane  is  essential, 
which  is  now  in  progress. 


5.  Conclusions 

p(  O2)  at  the  GDL  and  the  upper  channel  wall  at  the  cathode  of  an 
operating  PEFC  was  visualized,  revealing  the  different  distributions 
at  different  locations  in  the  gas-flow  channel.  At  low  Uq2,  p(02)  on 
the  upper  channel  wall  gradually  decreased  along  the  length  of  the 
gas-flow  channel,  whereas  that  at  the  GDL  was  found  to  decrease 


(a) 


Uo2  =  30%  (0.22  A  cnr2) 


(b) 


60%  (0.43  A  cm-2) 


(c) 


80%  (0.58  A  cm-2) 


Fig.  11.  Numerically-simulated  distributions  of  p(02)  on  the  cathode  catalyst  layer  (a)-(c),  the  local  current  density  on  the  MEA  (d)-(f),  and  the  water  vapor  activity  in  the  polymer 
electrolyte  membrane  (g)-(i).  The  cell  operation  conditions  used  for  the  simulation  are  the  same  as  those  shown  in  the  captions  of  Fig.  3.  (For  interpretation  of  the  references  to 
color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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significantly  at  the  center  of  the  length  of  the  flow  channel.  By 
increasing  U02 ,  the  steep  decrease  in  p( O2)  at  the  GDL  shifted  closer 
to  the  inlet.  The  change  in  the  distribution  of  p(02)  should  be 
reflecting  the  reaction  distribution  in  the  MEA,  which  was  sup¬ 
ported  by  the  numerical  simulation.  This  dual  visualization  tech¬ 
nique  combined  with  the  simulation  is  expected  to  be  used  for  the 
cell  and  MEA  designing. 
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List  of  Symbols 

i:  local  current  density  [A  m-2] 
io:  exchange  current  density  [A  m-2] 

6:  rate  of  catalyst  utilization  [-] 

G:  energy  of  activation  [J  kg-1] 

F;  Faraday  constant  [C  mol-1] 

R:  gas  constant  [J  mol-1  K-1] 

aa:  electrochemical  kinetic  rate  [-] 

ri a:  activation  over-potential  [V] 

r\x:  Ohmic  over-potential  in  electrolyte  membrane  [V] 

T:  temperature  [K] 

p:  density  [kg  m-3] 

h:  enthalpy  [J  kg-1] 

u:  velocity  [ms-1] 

p;  pressure  [Pa] 

s:  water  saturation  [-] 

pc;  capillary  pressure  in  porous  media  [Pa] 

av:  water  vapor  activity  [-] 

K:  permeability  [m-2] 

p:  viscosity  [Pa  s] 

r.  porosity  of  porous  media  [-] 

w;  molecular  weight  [kg  mol-1] 

c:  molecular  concentration  [mol  m-3] 

D;  gas  diffusivity  [m2  s-1] 
k:  thermal  conductivity  [W  m-1  K-1] 

mp:  mass  production  rate  by  electrochemical  reaction  [kg  m-3  s-1] 
rhgl:  mass  production  rate  by  phase  change  [kg  m-3s-1] 
c L>gi :  phase  change  rate  [s-1] 

de:  equivalent  hydraulic  diameter  of  gas-flow  channel  [m] 
nd'.  osmotic  drag  coefficient  [— ] 

Dw:  water  diffusivity  [m2  s-1] 

EW:  equivalent  weight  [kg  mol-1] 

A;  water  content  [— ] 

Superscripts 

ref:  reference  properties 
eff:  effective  properties 
k:  gas  component 

Subscripts 

(k):  =(a)  anode  side/=(c)  cathode  side  in  the  catalyst  layer 
(m):  polymer  electrolyte  membrane 
g:  gas  phase 
l:  liquid  phase 
s:  solid  phase 


